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Analyse van een object-dataset op een meer-dimensionale ru irate. 



De uitvinding heeft betrekking op een werkwijze voor het analyseren van een 
object-dataset op een meer-dimensionale ruimte. 



5 Zo 'n object-dataset representeeit 66n of raeer eigenschappen van het te 

onderz jeken object. Met name heeft de objecl-dataset betrekking op de dichtheidsverdeling in 
het te onderzoeken object en dan zijn de datawaarden de locale dichtheidswaaiden van (een 
deel van) het te onderzoeken object. De datawaarden kunnen bijvoorbeeld ook betrekking 
hebber op de verdeling van de temperatuur of de magnetisatie in het object. De rneer- 

10 dimemionale ruimte is raeeatal de drie-dimenstonale ruimte. Dan hebben de datawaarden 
betreklang op een volume verdeling van de betreffende eigenschap, bijvoorbeeld de 
dichtheiidsverdeling in een volume van het Le onderzoeken object De meer-dimensionale 
ruimte kan ook twee-dimensicmaal zijn. In dit geval hebben de datawaarden betrekking op een 
verdeling van de betreffende eigenschap in een vlak door het object, bijvoorbeeld de 

1 S dichtheuds verdeling Langs een dwarsdoorsnede door het object. 

De object-dataset kan op allerlei raanieten worden verkregen. In het bijzonder 
heeft de object-data set betrekking op een te onderzoeken patient. Zo'n objecc-dataset kan 
worden opgenomen met diverse technieken, zoals met 3D rOntgen rotationele angiografie, met 
computer tomografie of met magnetische resonantie afbeelding, of met magnetischc resonantie 

20 angiografie. 



Het artikel '3D rotational angiography: Clinical value in endovascular 
treatm*mt* in Medica Mundi 42(1998) door J. Moret et al. gaat over het bepalen van de 
25 doorsnsde van een bloedvat. 



Volgens de bekende werkwijze wordt uit de object-dataset een min of meer 
buisvorrnige structuur geidentificeerd die betrekking heeft op bloedvat van de te onderzoeken 
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patient, Deze ldentificatie gebeurt aan de hand van verschillen tussen datawaarden die 
betrekldng hebben op het bloedvat en op omringend weefsel. Dergelijke verschillen worden 
b reiki door bij het opnemen van de object-dataset aan de te onderzoeken patitat een geschikt 
contrar.tmiddel in het vaatstelsel toe te dienen. Volgens de bekende werkwijze wordt op het 
oog een snedeoppervlak min of meer dwars op de buisvorraige structuur in de object-dataset 
aangelisgd, Vervolgens wordt uit de dwarsdoorsnede van de buisvonnige structuur langs het 
snedeoppervlak de diameter van het betreffende bloedvat berekend. 

Een probleem bij de bekende werkwijze is dat het lastig is om het 
snedeoppervlak nauwkeurig loodrecht op de buisvonnige structuur te orienteren. Deze 
onnauvkeurigheid heeft ernstige fouten in de berekende diameter van het betreffende bloedvat 
tot gevDlg. 

Een doel van de uitvinding is een werkwijze te verschaffen waarmee analyse 
van eeti buisvonnige structuur in de object-data nauwkeuriger result aten oplevert dan de 
bekende werkwijze* Met name is een doel van de uitvinding een werkwijze te verschaffen 
waarmee de diameter van zo'n buisvonnige structuur nauwkeurig kan worden bepaald, 

Dit doel wordt volgens de uitvinding bereikt met een werkwijze voor het 
analyseren van een object-dataset waarin een buisvonnige structuur optreedt, 

- welke object-dataset datawaarden toevoegt aan posities in een meer-dimensionale ruimte, 
welke datawaarden beetrekking hebben op een te onderzoeken object, waarbij 

- een beginplaats in op of nabij de buisvonnige structuur wordt gekozen 

- gradiSnten op het oppervlak van de buisvonnige structuur en ter plaatse van de beginplaats 
woxden afgeleid 

- een snedeoppervlak wordt afgeleid door de beginpositie met richting zoveel mogelijk 
eve:iwijdig aan de gradignten in posities in een omgeving van de beginplaats en op het 
oppervlak van de buisvonnige structuur en 

- een dwarsdoorsnede van de buisvormige structuur langs het snedeoppervlak wordt 
afgeleid 

De object-dataset kan bijvoorbeeld worden verkregen met verschillende 
technielxn zoals magnetische resonantie angiografie, computer tomografle of 3D rotationele 
rBntgen angiografie. Met dergelijkc technieken wordt een object-dataset met datawaarden 
verkreg*n die de structuur van (een deel van) het vaatstelsel van de te onderzoeken patifini 
vertegeuwoordigt. 
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De buisvormige structuur betreft ecn gedeelte van de object-dataset waar de 
datawaarden afwijken van de omgeving van de buisvormige structuur. De afwijkende 
data* aarden liggen bijvoorbeeld in ecn geselecteerd bereik van datawaarden die groter zijn 
dan •» vooraf bepaalde drempel. of datawaarden die kleiner zijn dan een vooraf bepaald 
plafond. De posities in de meer-dimensionale ruimte waarbij de daiawaarden zo'n afwykende 
waste e hebben liggen langs een langgerekt min of meet buisvormig gebied. Dit buisvormig 
gebied kan de vorm hebben van een lege holla of een gevulde buis. De buisvormige structuur 
kan met segmentatietechnieken in de object-dataset worden aangeduid waarbij de posities met 
de afvijkende datawaarden worden geifdentificeerd. 

De beginplaats duidt de plaats aan waar de dwarsdoorsnede van de buisvormige 
structuur gewenst wordt. De beginplaats kan bu voorbeeld aangeduid worden door die in een 
weergave van de object-dataset op een beeldscherm aan te wijzen. 

De gradienten betreffen locale verschillen tussen datawaarden ter plaatse van 
het oppervlak. Volgens de uitvinding wordt de normaalvector op het snedeoppervlak 
nauwkeurig dwars op de gradienten op het oppervlak gericht. De normaalvector wordt zo 
gericlr dat de normaalvector binnen een nauwe tolerantie loodrecht op de gradienten staat. Zo 
lukt het dat het snedeoppervlak nauwkeurig nagenoeg loodrecht op buisvormige structuur 
wordt georienteerd. Merdoor is het mogelijk om uit de dwarsdoorsnede van de buisvormige 
structuur nauwkeurig de diameter van de buisvormige structuur af te leiden. Deze diameter is 
bijvoojbeeld een nuttig technisch resultaat bij het ondcrzoek van het stelsel van bloedvaten 
van de te onderzoeken patient. Met name lukt het om vemauwingen (stenoses) of uitatulpingen 
(an urismen) nauwkeurig te onderzoeken. 

Deze en andere aspecten van de uitvinding worden verder uitgewerkt aan de 
hand v«n de volgende uitvoeringsvoorbeelden die zijn gedefinieerd in de afhankelijke 
25 Conclusies. 

De normaalvector op het snedeoppervlak wordt bijvoorbeeld gericht door het de 
gewog«n som van de kwadraten van de inprodukten (scalair product) van de normaalvector en 
de gradifinten op het oppervlak van de buisvormige structuur in een (Weine) omgeving van de 
beginplaats te minimaliseren bij variatie van de normaalvector op het snedeoppervlak. De 
30 gewichien in de gewogen som kunnen bijvoorbeeld worden gebruikt om naar de invloed van 
gradient in posities dichterbij de beginplaats op de normaalvector op het snedeoppervlak 
groter t<; maken dan de invloed van gradienten in posities die wat verder van de beginposide 
liggen. 



# 

4 01.10.1999 



Bij voorkeur wordt uit posities in het snedeoppervlak met datawaarden in het 
gesel««teerde bereik een locaal centrum van de buisvormige structuur in het snedeoppervlak 
afgekid. Dit locale centrum vertegenwoordigt het midden van de buisvormige structuur in de 
dwarsdoorsnede langs het snedeoppervlak. Bijvoorbeeld Hgt het locale centrum op de positie 
5 die ongeveer gelijke afstanden heeft tot vrijwel alle posities in het snedeoppervlak met 
datawaarden in het geselecteerde bereik. Het blijkt dat in veel gevallen verreweg de meeste 
van dtae posities met datawaarden in het geselecteerde bereik en in het snedeoppervlak 
plaau>;n op de rand van de buisvormige structuur vertegenwoordigen. Bij voorkeur worden 
voor het afleiden van het locale centrum posities gebruikt waarin de gradienten magnitude 
10 hebben die groter zijn dan een minimumwaardc. Hierdoor wordt een nauwkeurige selectie van 
posities op de rand van de buisvormige structuur uitgevoerd. 

Een nog nauwkeuriger resultaat voor het locale centrum wordt op iteratieve 
wijze l>ereikt volgens de werkwijze zoals gedefinieerd in Conclusie 4. Bijvoorkeur wordt een 
verdelmg in viw gelijke grote sectoren gebruikt. Door per sector de minimum afstand van 
15 posities met datawaarden in het geselecteerde bereikt te bepalen worden posities buiten de 
bet reff snde buisvormige structuur uitgesloten. Wanneer in verschillende sectoren 
versch llende minimum afstanden tot de actuele centrumpositie van locale centrum optreden 
ligt de actuele centrumpositie niet nauwkeurig in het midden van de buisvormige structuur. 
Door de centrumpositie iets te vcrschuiven in de richting van de sector waar de grootstc 
20 minimum afstand optreedt wordt een nauwkeuriger schatting van de centrumpositie verkregen. 
Iiefst wordt de centrumpositie verschoven over een afstand die overeenkomt met de helft van 
het verschil tussen de minimum afstanden in tegenover elkaar gelegen sectoren. 

Volgens de uitvinding is het ook mogelijk om de lengteaa van de buisvormige 
structuur te volgen. paartoe wordt een vervolgpositie afgeleid die verschoven is in de richting 
25 van de tiormaal op het snedeoppervlak in de beginpositie. Dan wordt een volgend 

snedcoiJpervlak afgeleid door de vervolgpositie waarbij te werk gegaan wordt volgens 
Conclude 1. In het volgend snedeoppervlak wordt opnieuw een locaal centrum van de 
dwarsdoorsnede door de buisvormige structuur langs het volgende snedeoppervlak. Door 
herhaaldtoepassen worden opeen volgende locale centra langs de buisvormige structuur 
30 afgeleicL De magnitude van de onderlinge verschuiving langs de buisvormige structuur is door 
de gebrjiker te kiezen. Naar gelang er kleinere verschuivingen bij het afleiden van steeds weer 
een voljgende vervolgpositie worden toegepast wordt de buisvormige structuur preciezer 
gevolgc . Bij voorkeur wordt voor de zo verkregen reeks dwarsdoorsneden door de 
buisvoruige structuur langs respectievclijke snedcoppervlakken een reeks diameters berekend. 
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Dit lcvert een goed inzicht in het verloop van de diameter van de buisvormige structuur langs 
dc lergteas van die structuur. Dit resultaat is een nuttf g hulpmiddel om de mate van stcnose in 
het vi.alstelsel van de tc onderzoeken patient te onderzoeken. 

Bij voorkeur wordt de uitvinding toegepast om een buisvormige structuur die 
5 een reservoir met daarop aangesloten toe- en afvoerkanalen omvat. Het reservoir heeft typisch 
een aanzienlijk grotere diameter dan de toe- en afvoerkanalen. Zo'n siruatie doet zicfa 
bijvoorbeeld voor bij een aneurisma. Het eigenlijke aneurisma komt overeen met het reservoir, 
de toe- en afvoerkanalen komen overeen met de bloedvaten die Moed toe- en afvoeren naar en 
van hot eigenlijke aneurisma. Bij voorkeur wordt in een dergelijk situatie de werkwijze 
10 volgende Conclude 6 toegepast. Volgens de uitvinding lukt het om de toevoer- en 

afvoerpunten te bepalen door de toevoer- en afvoerkanalen te volgen volgens Conclusie 5. 
Liefst wordt dit volgen van de toevoer- en afvoerkanalen door de gebruiker op een monitor 
gevoljjd en kan de gebroiker het volgen becindigen zodra de aansluitingen op het reservoir zijn 
bereikt. Vervolgens wordt een middenpositie in het reservoir bepaald. Bijvoorbeeld wordt een 
15 deel van dc object-dataset gekozen welke het reservoir omvat en ten opzichte van het reservoir 
een klsin aantal posities buiten het reservoir. Het blijkt dat het zwaartepum van de 
dataw jarden in het geselecteerde deel van de object-dataset goed met het midden van het 
r servsir overcenkomt. Het snedeoppervlak door de toevoer- en afvoerpunten en door het 
midden van het reservoir scheidt het reservoir van de toevoer- cn afvoerkanalen. Het volume 
20 van he t reservoir aan de zijde van het snedeoppervlak dat naar de middenpositie is gericht 
gecft xiauwkeurig het effectieve volume van het reservoir weer. In het geval dat het reservoir 
met de toevoer- en afvoerkanalen een aneurisma met de daarop aangesloten bloedvaten 
vertegsnwoordigt geeft het effectieve volume van het reservoir een nauwkeurig resultaat voor 
de groDtte van het aneurisma hetgeen van belang is voor de verdere behandeling van het 
25 aneuriuna. 

Bij voorkeur wordt de werkwijze uitgevoerd met een daartoe geschikt 
geproixammeerd werkstation. Zo'n werkstation is bijvoorbeeld geladen met een 
computerprogramma met instructie voor het verrichten van de stappen zoals gedefinicerd in de 
Conclusies 1 tot en met 7. 

Een andere toepassing van de uitvinding op een buisvormige structuur met een 
reservoir met daarop aangesloten toevoer- en afvoerkanalen is gedefinieerd in Conclusie 7. 
Ook hierbij worden de toevoer- en afvoerpunten bepaald door de toevoer- en afvoerkanalen te 
volgen volgens Conclusie 5. In deze toepassing van de uitvinding wordt eon verbindingakanaal 
geinterpoleerd tussen het toevoer- en afvoerpunt door het reservoir heen. Vervolgens wordt het 
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verbindingskanaal van het reservoir afgezondcrd. Het resterende residu blijkt ecn goedc 
benadering van het effectieve reservoir dat goed bruikbaar is om de grootte van het aneurisma 
te bepslen. 

Deze en andere aspecten van de uitvinding warden bij wijze van voorbeeld 
5 nader toegelicht aan de hand van de volgende uitvoeringsvoaibeelden en aan de hand van de 
bijgevoegde tekening waarin 



Figuur 1 een voorbeeld toont van een object-dataset met daarin een buisvormige 
10 structuur waarin de toepassing van de uitvinding wordt aangeduid en 

Figuur 2 een voorbeeld van een stroomschema van de werkwijze volgens de 
uitvinding toont. 



15 Figuur 1 toont een voorbeeld van een object-dataset met daarin een buisvormige 

structuur waarin de toepassing van de uitvinding wordt aangeduid. Figuur 2 toont een 
voorbeeld van een stroomschema van de werkwijze volgens de uitvinding. Met name toont 
Figuur 1 de buisvormige structuur 1, in dit geval een weergave van een deel van het vaatstelsel 
van de te onderzoeken patient, die al uit de eventucel groterc object-dataset is gesegmenteerd. 

20 Dit is iiangeduid als de stap 10 in Figuur 2. De segmentate komt er bijvoorbecld op neer dat 
alle daawaarden op een vaste waarden worden gezet in posities waarin de oorspronkehjke 
object-dataset datawaarden buiten het geselecteerd bereik heeft. Het vaatstelsel omvat een 
aneurisma 3, dat als reservoir fungeert, ecn bloedvat 2 dat bloed naar het aneurisma 3 toevoert 
en een bloedvat 4 waarmee bloed uit het aneurisma 3 wordt afgevoerd. Dus de bloedvaten 2 en 

25 4 fung-aen als toevoer- en afvoerkanaal. Het beginpunt (BP) wordt door de gebruiker gekozen, 
bijvoorbeeld door het met behulp van een buis en/of een toctsenbord met een cursor op het 
beeldscherm van een werkstation aan de wijzen. Om het beginpunt (BP) wordt een bol (Bl) 
genomen zoals aangeduid in stap 20. Op punten op de rand van de buisvormige structuur in de 
bol (B I) worden de gradienten (gij) van de datawaarden berekend, dit is aangeduid in stap 30. 

30 Vervolgens wordt in stap 40 met een minimalisatie proces de normaal vector n, berekend, zo 

dat de jom J • « u )* mi™maal is. De gewichtsfactoren w } zijn door de gebruiker in te 

stellen liefst nemen de gewichtsfactoren af naar gelang de afstand van de positie waar de 
gradient glj berekend wordt tot het beginpunt BP groter is. Het snedeoppervlak SN1 heeft 
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norn^alvector /i, en gaat door het beginpunt BP. De dwarsdoorsnede door de buisvormige 
stmctuur 1 langs het snedeoppervlak SN1 wordt genoroen in stap 50. Uit deze dwarsdoorsnede 
wordt de locale diameter rfvan dc buisvormige structuur 1. d.w.z. van het bloedvat, afgeleid in 
stap 51. Verder wordt in stap 52 uit de dwarsdoorsnede de locale centrumposttie van de 
buisvc .rmige strucruur afgeleid door de iteratieve procedure volgens Conclusie 4. 

Om vanuit het beginpunt BP de buisvormige structuur te vervolgen wordt in 
stap 60 de vervolgpositie VP afgeleid in de richting van de normaalvector n,, die locaal 
nauwkeurig langs as van het bloedvat is gericht De grootte van de verschuiving wordt naar 
gelang te kronkeligheid van de buisvonnige sttuctuiff doorde gebwiker ingesteld. m Figuur 1 
is een vervolgpositie VP getoond die wordt bereikt na een groot aantal verschuivingen lungs 
de as van het bloedvat 2. Wanneer de vervolgpositie is nog niet de door de gebruiker gewenste 
eindpositie (EP) is worden de stappen 20.30,40,50,51.52 en 60 herhaald. Zo wordt langs het 
betreff sndc bloedvat de as bepaald en de diameter langs de as. 

Bijvoorbeeld is het eindpunt EP het toevoerpunt TP van het bloedvat 2 naar het 
aneuri* ma 3. Door vanaf een ander beginpunt BP' op bloedvat 4 het bloedvat 4 te vervolgen 
volgem de procedure die in figuur 2 is aangeduid wordt het afvoerpunt AP bereikt. De 
gebruil:er lean bijvoorbeeld zelf het vervolgen van het bloedvat afbreken wanneer naar het 
anatorr isch inzicht van de gebruiker het toevoer- of afvoerpunt is bereikt. Het is ook mogelijk 
het vervolgen van het bloedvat af te breken wanneer de locale diameter van het bloedvat 
plotseling sterk toeneemt. Door het toevoer- en afvoerpunt (TP, AP) wordt een verbindingslijn 
/ getrol Jccn. Verder wordt de middenpositie Z van het aneurisma 3 berekend als het 
^ z waarttspunt aan de hand van de datawaarden in het gebied O. Het gebied G is bij voorbeeJd 
een bol om het aneurisma 3. Het blijkt dat goede resultaten voor de middenpositie Z worden 
bereikt wanneer ervoor gezorgd wordt dat het gebied G niet al te veel punten in de 
25 gesegmenteerde dataset bevat buiten het aneurisma. De normaal vector m is gericht 

tegengesteld aan de loodlijn van uit Z op de verbindingslijn I. Vervolgens wordt nog een 
snedeojipervlak SN3 bepaald door de verbindingslijn en met normaalvector A . Het blijkt dat 
dit sneceoppervlak SN3 nauwkeurig de toevoer- en afvoerkanalen, in dit voorbeeld de 
bloedveten 2 en 3. van het aneurisma 3 scheidt. 
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1. Een werkwijze voor het analyseren van era objcct-dataset waarin een 
buisvormige structuur optreedt, 

- wclke object-dataset datawaarden locvoegt aan posities in een meer-dimensionale ruimte, 
wclke datawaarden beetrekking hebben op een te ooderzoeken object, waarbij 

5 - een beginplaats in op of nabij de buisvormige structuur wordt gekozen 

- gradienten op het oppervlak van de buisvormige structuur en ter plaatse van de beginplaats 
worden afgeleid 

- een snedeoppervlak wordt afgeieid door de beginpositie met richting zoveel mogelijk 
evenwijdig aan de gradi&nten in posities in een omgeving van de beginplaats en op het 

10 ojpervlak van de buisvormige structuur en 

- een dwarsdoorsncde van de buisvormige structuur langs het snedeoppervlak wordt 
af jeleid. 

2. Een werkwijze volgens Conclusie 1, waarin de richting van het snedeoppervlak 
IS wordt afgeleid door minimalisatie van de som van de kwadraten van de inprodukten van de 

normt alvector op het snedeoppervlak en de gradifetten. 

3. Een werkwijze volgens Conclusie 1, waarin een locaal centrum van de 
buisvormige structuur in het snedeoppervlak wordt afgeleid uit posities in het snedeoppervlak 

20 met djitawaarden in een vooraf gekozen bereik, waaibij zo'n bereik in het bijzonder 
datawaarden betreft groter dan een drempelwaarden of datawaarden kleiner dan een 
plafor dwaarde. 



4. Een werkwijze volgens Conclusie 3, waarbij het snedeoppervlak in meerdere 

25 sectors verdeeld wordt en 

- een schatting van een centrumpositie van het locale centrum van de dwarsdoorsnede van 
de buisvormige structuur in het snedeoppervlak wordt gemaakt 

- voor de respectievelijke sectoren een minimum afstand van posities met datawaarden in 
het vooraf gekozen bereik tot een actuele geschatte centrumpositie wordt afgeleid 
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- een maximum van de minimum afstanden bij de respectievelijke sectoral wordt bepaald 
~ een nieuwe schatting van de centrumpositie wordt afgeleid, warbij dc nieuwe schatting van 

de centrumpositie ten opzichte van de actuele schatting is verachoven in het 
snedeoppervlak in de richting naar de sector waarin het maximum van de minimum 
5 af itanden ligt en over een afstand die afhangt van de vcrschillen tussen posities waarbij in 
tejjenover elkaar gelcgen sectoren de minimum afstanden optreden 

- eventueel de stappen b,c en d worden herhaald waarbij de nieuwe schatting van de 
centrumpositie als actuele schatting wordt gebruikt. 

10 -5: 



Een-werkwijzevdlgens Conclusie 3, waarbij 

- een vervolgpositie wordt bepaald welke verschoven is ten opzichte van de beginpositie in 
de richting van de normaal op het snedeoppervlak door de beginpositie en 

- een volgend snedeoppervlak door de vervolgposide wordt afgeleid, waarbij de richting van 
dii volgend snedeoppervlak zoveel mogelijk cvenwijdig is aan de gradi*nten in eft volgend 

15 sn ideopperv] ak en op het oppervlak van de buis vonnige structuur. 

6 Een werkwijzc voor het analyseren van een object-dataset waarin een 

buisvormige structuur optreedt, 

- w<lke object-dataset datawaarden toevoegt aan posities in een meer-dimensionale ruimte, 
20 wj.arbij 

- in de buisvormige structuur een tocvoerkanaal. een afvoerkanaal en een reservoir optreden, 
en de toevoer- en afvoerkanalen op het reservoir zijn aangesloten en 

- een aanvoerpunt in het toevoerkanaal ter plaatse van de aansluiting van het toevoerkanaal 
op het reservoir wordt bepaald, 

25 - ecu afvoerpunt in het afvoerkanaal ter plaatse van de aansluiting van het afvoerkanaal op 
hei reservoir wordt bepaald, 

- ecu middenpositie van het reservoir wordt bepaald en 

- een snedeoppervlak wordt afgeleid door de aan- en afvoerpunten, waarbij de normaal op 
hei snedeoppervlak langs de loodlijn uit de middenpositie op de lijn door de aan- en 

30 afvoerpunten staat. 

7 - Een werkwijze voor het analyseren van een object-dataset waarin een 

buisvormige m^amaptneSt, " ~ 
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- we Ike object-dataset datawaarden tocvoegt a an posities in een meer-dimensionale ruimte, 
waarbij 

~ in de buisvormige structuur een toevoerkanaal, een afvocrkanaa] en een reservoir optreden, 
en de toevoer- en afvoerkanalen op het reservoir zijn aangesloten cn 
S - een aanvoerpimt in het toevoerkanaal ter plaaise van de aansluiting van het toevoerkanaal 
op het reservoir wordt bepaald, 

- een afvoerpunt in het afvoerkanaaJ ter plaatse van de aansluiting van het afvoerkanaal op 
he1 reservoir wordt bepaald, 

- een verbindingskanaal wordt gelnterpoleerd tu&sen de aan* en afvoerpunten door het 
10 reservoir en — 

* een residu wordt afgeleid als verschil tussen het reservoir en het geYnteipoleecde 
verbindingskanaal. 
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The method relates to analysis of a tube -like structure 1, such as a rendering of 
a patient's vascular system. In a begin point (BP) gradients g }j to the surface of the tube-like 
struetire 1 are computed Next a normal n essentially perpendicular to the gradients is derived 
and a cross section is taken along a cutting plane having n as its normal. A subsequent 
5 position (VP) is taken by performing a small shift in the direction of h . By repetition of the 
procedure the axis of the vessel is tracked. 



Figure 1 
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